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We investigate scaling of work and efficiency of a photonic Carnot engine with the number of 
quantum coherent resources. Specifically, we consider a generalization of the “phaseonium fuel” for 
the photonic Carnot engine, which was first introduced as a three-level atom with two lower states 
in a quantum coherent superposition by [M. O. Scully, M. Suhail Zubairy, G. S. Agarwal, and H. 
Walther, Science 299 , 862 (2003)], to the case of -|- 1 level atoms with N coherent lower levels. 
We take into account atomic relaxation and dephasing as well as the cavity loss and derive a coarse 
grained master equation to evaluate the work and efficiency, analytically. Analytical results are ver¬ 
ified by microscopic numerical examination of the thermalization dynamics. We find that efficiency 
and work scale quadratically with the number of quantum coherent levels. Quantum coherence 
boost to the specific energy (work output per unit mass of the resource) is a profound fundamen¬ 
tal difference of quantum fuel from classical resources. We consider typical modern resonator set 
ups and conclude that multilevel phaseonium fuel can be utilized to overcome the decoherence in 
available systems. Preparation of the atomic coherences and the associated cost of coherence are 
analyzed and the engine operation within the bounds of the second law is verified. Our results bring 
the photonic Carnot engines much closer to the capabilities of current resonator technologies. 

PACS numbers: 42.50.Ar,05.70.-a,07.20.Pe 


I. INTRODUCTION 


A practical figure of merit to compare fuel and bat¬ 
tery materials is the specific energy, or energy to mass 
ratio EH- As a material constant, it measures the en¬ 
ergy that will be harvested by using a unit mass of the 
material. About a decade ago, a highly non-traditional 
fuel, called “phaseonium”, which is a three level atom 
with two lower states in a quantum coherent superpo¬ 
sition, was proposed to be used in a Photonic Carnot 
engine (PCE) [^. Phaseonium engine could work with a 
single heat bath and a phaseonium reservoir This 

proposal stimulated much interest to quantum heat en¬ 
gines [sj-Eoj]. It was later argued that existing resonator 
systems can not implement such an engine, due to high 
cavity losses and atomic dephasing [^. Here, we ad¬ 
dress a fundamental question of how the specific energy 
of phaseonium fuel is scaled with the number of quantum 
coherent levels. A favorable scaling law against decoher¬ 
ence and dephasing could bring the phaseonium engine 
closer to available practical systems. 

We describe multilevel generalization of phaseonium 
fuel in Fig. [TJ The block diagonal density matrix p of 
an Af -b 1 level atom is shown in Fig. [IJa). The excited 
level, denoted by “a”, and the lower levels, denoted by 
“6i, 62,.., &Ar” are well separated from each other by an 
energy Q measured from the central lower level 6Ar/2 as 
shown in Fig. [IJb). The lower levels can be degenerate 
or non-degenerate. The diagonal elements paa, and 
with b G {61,62,.., &Ar}, determine the level populations, 
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while the off diagonal elements pbb', with b' ^ b, indi¬ 
cate the coherence between the levels. Coherence can be 
characterized by the magnitude and phase of the complex 
number pbv ■ 

Though both the amplitude and the phase of coher¬ 
ent superposition states can be controlled in experiments 
(2^, the main control variable for the photonic Carnot 
engine is the phase of the coherence as the amplitude 
is required to be small enough to keep the system only 
slightly out of thermal equilibrium. The complete graphs 
in Fig. [Ija) have N nodes and N{N — l)/2 links, repre¬ 
senting the atomic energy levels and the coherences be¬ 
tween them, respectively. The simplest graph has N = 2 
nodes, which is the case of the original phaseonium pro¬ 
posal . The interplay between quantum coherence and 
energy discussed in photon Carnot engine revealed 
that the energy content of the phaseonium with N = 2 
can be optimized at a certain phase of the coherence. 
We could envision as if we are considering more com¬ 
plex, larger, phaseonium molecules with the graphs hav¬ 
ing N > 2, corresponding to A^-bl level atom phaseonium 
(NLAP). 

We can imagine different phaseonium molecules can be 
possible for a given atom of unit mass and explore how 
the specific energy of the atom depends on the size of the 
phaseonium molecule characterized by N. Next to the 
phase of coherence, N becomes another control param¬ 
eter which could favourably contribute to the enhance¬ 
ment of the specific energy of the single atom quantum 
fuel. For ^ 1, the number of coherences would scale 
quadratically, ^ N^. If the quadratic coherence scaling 
could be translated into the energy content of the atomic 
fuel, we could overcome the cavity losses for implementa¬ 
tion and boost the performance of quantum Carnot en- 
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Figure 1. (Color online) N+1 level atom phaseonium (NLAP) 
fuel, a) Density matrix p and complete graph representations 
of NLAP. p is + 1 dimensional square matrix. Its coherent 
block can be represented by a complete graph with N nodes 
and S = N{N — l)/2 links. Graphs were shown up to A = 5 
number of nodes, b) NLAP for non-degenerate and degener¬ 
ate atoms. The excited state is denoted by a and the lower 
levels are denoted by bi with i = 1..N. The upper level is well 
separated from the lower levels by an energy Ml measured 
from the central lower level 6 iv/ 2 - 
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Figure 2. (Color online) Photonic Carnot engine with A -|- 1 
level atom phaseonium (NLAP) fuel. Photon gas in a high 
quality cavity of frequency is the working substance and 
the mirrors of the cavity play the role of the piston. NLAP 
leaves the hohlraum at temperature Th and is subsequently 
prepared in a state with quantum coherence among its lower 
levels characterized by A(A—1)/2 phase parameters (f>ij with 
i,j = Created NLAPs are repeatedly injected into the 

cavity at a rate r in the quantum isothermal expansion pro¬ 
cess, where heat Qin is transferred to the cavity. The cycle 
continues with quantum adiabatic expansion and quantum 
isothermal compression and is completed with a quantum adi¬ 
abatic compression. An amount of heat Qout is rejected into 
the entropy sink in the isothermal compression. 


gine for applications. From fundamental point of view, 
such scaling analysis could reveal profound difference of 
quantum fuel from a classical resource as such a scal¬ 
ing cannot exist without quantum coherence. Complete 
graphs of phaseonium molecules serve more than a simple 
counting of coherences. They emphasize the generality 
of our question we address in the present contribution. 
Can we beat decoherence with the scaling advantage of 
quantum coherent resources? 

We have recently numerically examined a superradi¬ 
ant quantum Otto engine which exhibits a similar 
scaling law of work with the number of atoms; though 
the efficiency is independent of the temperature and the 
number of atoms. Our present results are fully analyti¬ 
cal, and for a Carnot cycle which can harvest work from 
a single heat bath and the scaling laws with the number 
of coherences are both for the work and the efficiency. 
Furthermore, specific energy of a superradiant cluster in¬ 
crease linearly with the number of atoms, while here for 
the single atom phaseonium fuel it is quadratic. Quan¬ 
tum advantage in the charging power of quantum bat¬ 
teries with the number of qubits has been examined very 
recently 23|. Linear scaling of work with the number of 
qubits is reported; while due to a quantum speed up of 
the ODCTation time, the charging power scales quadrati- 
cally 23| . The preparation of phaseonium fuel and using 


it in PCE can be compared to charging and discharg¬ 
ing a single qubit quantum battery with multiple quan¬ 
tum coherences. Phaseonium fuel or battery allows for 
quadratic scaling in harvested work, efficiency, and the 
specific energy with the number of quantum coherences. 
We examined the charging or preparation cost of the 
phaseonium battery and compared it with the harvested 
work by PCE. Our result verified that the second law is 
obeyed in our system. 


II. NLAP MODEL AND SYSTEM DYNAMICS 

The operation of photonic Carnot engine is described 
in Fig. [5] The working fluid of the engine is the photon 
gas in a high quality single mode cavity of frequency LI. 
The radiation pressure by the cavity photons applies on 
the cavity mirrors playing the role of the piston of the 
engine. The quantum fuel of the engine is an NLAP. The 
quantum Carnot cycle consists of two quantum isother¬ 
mal and two quantum adiabatic processes. 

In the isothermal expansion, NLAPs are generated and 
injected into the cavity at a rate r. The interaction 
time T between an NLAP and the cavity field is short, 
T < 1/r, so that only one NLAP can be present in the 
cavity [^. Coherences in NLAP are characterized by 
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N{N — l)/2 phase parameters with i,j = 1, 2, N. 
Coherent superposition states in iV +1 level atom system 
can be gene rated by stimulated Raman adiabatic pas¬ 
sage , Morris-Shore transformation 12711 , or quan¬ 

tum Householder reflection techniques [2^,[2j Thermal- 
ization of the single atom can be considered relatively fast 
and hence the injection rate would be limited by the time 
of coherence preparation. The choice of specific technique 
of coherence induction depends on the details of a partic¬ 
ular implementation. If the amplitudes of the coherences 
are much smaller than the level populations, then NLAP 
can be assumed in an approximate thermal equilibrium 
with a thermal reservoir (hohlraum) at a temperature T^. 
During the interaction, the mean number of photons, h, 
and the cavity temperature increases; while the expan¬ 
sion cools down the cavity when there is no atom inside. 
Repeated injection of NLAPs into the cavity maintains 
the cavity field at a temperature by transferring a 
total amount of heat into the cavity as Qin. is an ef¬ 
fective temperature defined in terms of the steady state 
photon number fitf, as = fiD/fcln(l -f l/h^), with k is 
the Boltzman constant. It can be higher than Th in the 
presence of coherence [^. The cavity volume, and hence 
the frequency, change negligibly, AD <C D. 

The cycle continues with an adiabatic expansion where 
the entropy remains constant and the temperature drops 
as the D changes appreciably. Following step is the 
isothermal compression in which heat Qout is transferred 
from cavity to a cold reservoir at a temperature Tc- The 
cycle is completed by adiabatic compression where the 
temperature is raised back to T^. 

The net work extracted from the cycle is IFnet = Qin — 
Qout where Qin — ^'Ud f^out — ^c(' 5^3 ‘ 5 ' 4 )' 

The mean photon number hi and the temperature Ti 
at the beginning of the stage determine the entropy 

51 by Si = fcln(hi -I- 1) -I- Mlfii/Ti. Using Si = 54, 

5 2 = 53 , Ti = T 2 = T^, and T 3 = T 4 = Tc, we write 

kFnet = (70 — Tc )(52 — 5i). The efficiency of the engine 
is defined as ry = Wnet/Qin- It reduces to 77 = 1 — Th/T^. 
This coincides with the standard definition of thermody¬ 
namic efficiency in Carnot cycle and used in the origi¬ 
nal proposal of the phaseonium fuel as well as in 

the arguments against its feasibility in the presence of 
decoherence channels [2l|. In order to present results 
comparable to the previous works, we calculate the effi¬ 
ciency as defined by these studies. To avoid any mislead¬ 
ing impressions however, we emphasize that in practical 
considerations round trip efficiency can be more relevant 
figure of merit. The round trip efficiency of the engine 
should include the cost of the preparation of the quan¬ 
tum coherent atom; which would ensure the validity of 
the second law Q . On the other hand, it was noted that 
the cost of quantum fuel can be expensive but it is 
still appealing as it can be used to harvest work from a 
single heat reservoir. Our objective here is not to discuss 
if such PCEs can be efficient enough for certain awlica- 
tions, but to examine if such devices, proposed in [^, can 
produce positive work in the presence of decoherence by 


exploiting a scaling advantage of multiple coherence re¬ 
sources, in contrast to the negative conclusions of earlier 
studies plj] . 

During the adiabatic process n does not change so that 
hi = h 4 = {exp{hn/kTc) — 1)~^ and h 2 = h^. These 
relations reveal that work and efficiency of the photonic 
Carnot engine can be calculated by determining the h^ 
at the end of the isothermal expansion stage. 

This manuscript organised as follws: In Sec. II, we de¬ 
scribe the N-level phaseonium model, and PCE system 
dynamics. We review and discuss the analytical and nu¬ 
merical verification of the analytical results in Sec. III. 
We also estimate the preparation cost of NLAP in this 
section. We conclude the results in Sec. IV. 


III. RESULTS AND DISCUSSIONS 

In order to find the h^, we solve h^ = ~ ® 

where pnn = Here \n) is the Fock number state 

for the cavity photons and p is the reduced density matrix 
of the cavity field. The equation of motion for p can 
be obtained by tracing the equation of motion of the 
complete system over atomic degrees of freedom 

= ( 1 ) 

k 

where = Hq -|- Hj is the Hamiltonian of the arbitrary 
k*^ atom in the interaction picture relative to the cavity 
photons, with Hq = tvjJa\a){a\ + 

Hj = |a)(6i|ae“®^‘ -I- H.c.. Here hua, hujbi are 

the energies of atomic states |a) and \bi), with i = 1..N, g 
is the coupling rate between the atom and the field, and 
a is the photon annihilation operator. The model Hamil¬ 
tonian describes a situation where N + 1 level atom is 
coupled to a single mode cavity in a fan shaped transition 
scheme. A more realistic model requires consideration of 
multimode cavity coupled to an atom with multiple up¬ 
per and lower hyperfine levels Bill. Such models can 
be reduced to effective single mode cavity and multilevel 
atom interactions [s^ or can be directly described by 
generalized master equations of micromasers [3^ . Atoms 
with fan shaped degenerate level schemes are also stud¬ 
ied from the perspective of generating large superposition 
states 1^. The central question for us here is the 
dependence of work and efficiency on the number of the 
superposed quantum states and we will only consider sin¬ 
gle upper level and a set of degenerate or non-degenerate 
lower levels for simplicity. 

Analytically calculating the right hand side of the 
Eq. ([T]), we find (see Appendix for details) 

Pnn = -rg‘^{KaPaa[in + l)pnn “ 

+ (Sfcl KbiPbibi + J2t<j \Pbibi I) 

X [np 

nn (n -I- l)p„+i_„+i]}. 


( 2 ) 
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where the coefficients Ka,Kb. and depend on the 
atomic relaxation rate 7 , atomic dephasing rate 7 ^, de¬ 
tuning parameter Ai = ujabi — with ujabi = bja — LOb^, 
and iObibj = — Wi, as well as the coherence parameters 

4>ij and I, by the relations given in the Appendix. 
Thus, we obtain the rate of change of average photon 
number 


'^(f) — ‘^9 {^aPaaij^cj) 1 ) (^) 

where Rg„ = X;*=i Kb.Pb.b, and i?g, = \pb,b^\. 

The equation of motion for the evolution of population 
elements in the density matrix coincides with the ther- 
malization dynamics of a resonator coupled to a heat 
bath. Accordingly, the coarse grained dynamics effec¬ 
tively describes sequence of NLAP injected into the res¬ 
onator as a mesoscopic ensemble of fV -|- 1 level atoms 
acting as a heat bath. The off diagonal elements of the 
density matrix or the coherences can be kept vanishingly 
small to describe the steady state approximately as a 
thermal equilibrium state. The corresponding effective 
temperature can be determined by the modified detailed 
balance condition to reach such a quasi equilibrium state 
in Eq. @ which gives 


Kg /-hn\ 

Rgo + Rgc \ ) 


( 4 ) 


The detailed balance between the thermal reservoir at Th 
and the photon gas in the resonator is broken but there is 
a modified detailed balance between the coherent atomic 
ensemble and the resonator photons. Accordingly the 
resonator can reach a thermal equilibrium at a different 
temperature than Th- 

The steady state of the Eq. ([3|) yields the average pho¬ 
ton number as 




n 


1 -I- n 


ifaPaa 


( 5 ) 


where h = {Rgg/Kgpgg — 1 )“^ is the average pho¬ 
ton number in the absence of coherence. Using = 
{exp{hQ/kT^) — 1 )“^, we determine the effective cavity 
temperature as 




Th 


1 I - fio 

1 -I- 




( 6 ) 


by using high temperature approximations fi^ Ri kT^/Ml 
and n r; kTh/Ml in Eq. ([5|). 

Therefore, the efficiency of the photonic Carnot engine 
in the case of NLAP becomes 


Tg_ R, 


'n<f> = 'nc- 


gc 


Th KaPa 


( 7 ) 


where rjc = ^ — X/Th is the Carnot efficiency. Note that 
for Tg = Th, r]c = 0 but rj^p could have a positive value for 
particular values of control parameters 4 )i, 4 > 2 , In 

order to get further analytical results we will make some 
simplifying assumptions. 


We focus on degenerate NLAP case to proceed analyt¬ 
ically, for which Eg = U, = 0,* = 1..N, ujgb^ = O, 
Ai = 0, oJbibj = 0 and Kg = 2N In addition, we 
consider phase locked equal amplitude coherences with 
(pij = (j) and \pbibj I = A. Hence the coefficients in 
Eq. ([2]) become = 4 cos(/)/ 7 ^, Rgg = 2NPg/j^ and 
Rg^ = 2N{N — 1) cos(/)A/ 77 , and hence Eq. ([3]) reduces 
to 


np, = 2pLN[{Pg - Pg + N^X)np, + Pg] - KUp, (8) 


for A > 1 , (^ = TT, where p = Pg = pgg = 

exp{-f3Eg)/Z, Pg = pb,bi = 1/^ with Z = exp(-,5ii;a) -I- 
N. Here we introduced k and with |^| < 1, as the de¬ 
coherence rate due to the dissipation in the cavity and a 
phenomenological decoherence factor due to atomic de¬ 
phasing, respectively [2l|. While the dephasing factor 
is phenomenologically introduced in 21 | we provide its 
rigorous microscopic derivation in the Appendix. 

Steady state solution of Eq. ([5]) yields an effective tem¬ 
perature given by Tp = Th/{1 + F{Th)) in the high tem¬ 
perature limit where 


= (9) 

with h = Pg/{Pg — Pg). Eor small coherence and deco¬ 
herence terms in F{Th), an approximate expression can 
be written for the effective temperature 

Tp = Th(^l + N^^Xn - 

This result shows that if the reduction of the magni¬ 
tude of coherence due to dephasing is slower than the 
quadratic increase with N, then the multilevel coher¬ 
ence could be used to beat the decoherence induced 
by the cavity dissipation. The magnitude of coherence 
A is limited by the positivity requirement of the den¬ 
sity matrix as well as the thermal equilibrium require¬ 
ment of the cavity field. The former condition require 
\Pbibj \ < {pbibiPbjbj)^^"^ so that A < 1/A for A > 1 as 
Pbibi 1/fV for A ^ 1. Accordingly one can fix the 
coherence amplitude A as a constant as long as it re¬ 
mains smaller than 1/A for the range of A values. Eor a 
realistic number of levels this is not a very restrictive con¬ 
dition. More severe limitation on A is due to the quasi¬ 
equilibrium condition of the photon gas. We will take 
A ^ 10“® and consider A < 40 in our numerical examina¬ 
tions. In the classical asymptotical limit of A —>■ oo then 
A —^ 0 as 1/A and hence the quadratic scaling reduces 
to a linear one for which the specific energy becomes a 
constant as expected for classical systems. Mesoscopic 
systems in quantum regime are therefore necessary to 
exploit the quadratic scaling in the specific energy. 

We note that the coarse grained dynamics is designed 
on purpose to determine the steady state by rapid conver¬ 
gence using a numerically efficient dynamical equations. 
Analytical solution of the mean photon number dynam¬ 
ics for the degenerate case and with A = 0 ,k = 0,^ = 1 


2p 


( 10 ) 











A. Analytical results for 
modern resonator systems 
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gives 


n^=n-{n-no)e ‘Z****, ( 11 ) 

where tth = — Pg) is the thermalization time. 

We will first discuss the predicted analytical state states 
by the coarse grained master equation in modern res¬ 
onator settings then examine the exact numerical desrip- 
tion of the dynamics of the system to verify the analytical 
results in the subsequent subsections. 



Figure 3. (Color online) Extracted work {W) and efficiency 
(rf) of photonic Carnot engine, with A -|- 1 level atom phaseo- 
nium (NLAP) fuel, depending on the number of degenerate 
coherent ground state levels N, for different decoherence fac¬ 
tor models (a)-(b) ^ = exp (—a;), (c)-(d) ^ = exp (—A®), (e)- 
(f) ^ = exp (—A^®), where a; = 7,^/7 . Coherence parameter 
is A = 10“® and the initial thermal coherent atomic temper¬ 
ature is Tji = 4 in units of ftfl/fcs. x values are 0.15, 0.1, 
0.05, 0.001 for (a)-(b), 0.14, 0.12, 0.1, 0.08 for (c)-(d) and 
0 . 012 , 0 . 01 , 0.008, 0.006 for (e)-(f) in decreasing order for the 
lower to upper curves respectively. The plots are given for 
the circuit QED parameters in [^ . The quantities g — 0.01, 
r = 1 X 10“"^, K = 6.25 X 10”“*, and 7 = 5 x 10“®, which are 
the coupling coefficient to the cavity field, atomic injection 
rate, cavity loss term, and atomic decay respectively, are di¬ 
mensionless and scaled with the resonance frequency fl 10 
GHz. 77 is dimensionless, W is dimensionless and scaled with 

n. 


In the high temperature limit (T 17), the en¬ 
tropy change in the isothermal expansion stage is AS' = 
fcAfl/n and the heat input becomes Qin = T/jAS. The 
work output at Th = is found to be IT = Qin?? where 
r] = — k/2^), respectively. In superconduct¬ 

ing circuit, microwave and optical resonators, it is es¬ 
timated that K/2fj,^X ^ 10 [^. should be replaced 
by A(A—1)/2 for smaller number of levels. Accordingly, 
by using five or more level quantum phaseonium fuel, the 
working flnid can beat quantum decoherence to harvest 
positive work. 

In Fig. [21 we plot the work output and efficiency 
of the photonic Carnot engine with degenerate NLAP 
fuel, depending on the number of quantum coherent lev¬ 
els. We consider N independent as well as N depen¬ 
dent scaling models [s^ for the decoherence factor and 
take ^ = exp (—a:) in Fig. |21a)-(b), ^ = exp(—Acc) in 
Fig.|^c)-(d), and ^ = exp(—A^a:) in Fig.|3De)-(f), where 
a; = 70/7 as shown in the Appendix. The plots are given 
for the circuit QED parameters in [2l|. We consider 
larger atomic dephasing rates than the typical values to 
demonstrate its limiting effect on W and rj. The plots in¬ 
dicate that even when there is large dephasing, which can 
increase with A linearly or quadratically, W and 7 can 
retain their quadratic power law with A up to a critical 
A. 

Similar results are found for the cases of optical and 
microwave cavities. We see both in Fig. |21a)-(b) and 
Fig. |5l(a)-(b) that when dephasing is independent of A, 
the work output and efficiency increases quadratically 
with the number of coherent levels. If the dephasing 
rate is increasing linearly with A as in Fig. |TKc)-(d) and 
Fig.[5Dc)-(d), or if it is increasing quadratically with A as 
in Fig.|3Ke)-(f) and Fig.|5Ke)-(f), the work output and ef¬ 
ficiency of the photonic engine is enhanced quadratically 
with the number of coherent levels only up to critical A. 
beyond which the work output and efficiency decays ex¬ 
ponentially due to the dominating effect of decoherence. 

B. Numerical verification of the theory 

In order to perform a faithful numerical simulation of a 
typical set up described in the theory, we investigate the 
injection process in detail. We assume a regular atomic 
injection of Rydberg atoms into a Fabry-Perot cavity [s^] 
with an atomic interaction time r with the cavity field 
and an empty cavity time tq such that 1 /r = t-\-tq where 
r is the injection rate. During the time interval r, the 
hamiltonian is 

N 

H = a;a|a)(a| -I- D,a^a + |a)(^i|a + H.c.) (12) 

i=l 

while for the time interval tq, H = (fi = 1 and 

Wf,. = 0 for degenerate ground state levels). 
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Figure 4. (Color online) Extracted work {W) and efficiency 
(rf) of microwave resonator photonic Carnot engine, with 
+ 1 level atom phaseonium (NLAP) fuel, depending on 
the number of degenerate coherent ground state levels N for 
different decoherence factors (a)-(b) ^ = exp(—x), (c)-(d) 
^ = exp(—A'"*), (e)-(f) ^ = exp(—A^^a;), where x = 7 ,^/ 7 . 
Coherence parameter is A = 10~® and the initial thermal co¬ 
herent atomic temperature is = 4 in units of hQ/ks- x 
values are 0.1, 0.15, 0.05, 0.001 for (a)-(b), 0.12, 0.1, 0.08, 0.06 
for (c), 0.14, 0.12, 0.1, 0.08 for (d), 0.12, 0.1, 0.08, 0.06 for (e) 
and 0.01, 0.008, 0.006, 0.004 in decreasing order for the lower 
to upper curves respectively. The parameters g = 9.21 x 10“^, 
r = 6.47 X 10"®, K = 1.96 x 10"® and 7 = 9.54 x 10”“ are the 
atom-field coupling coefficient, atomic injection rate, cavity 
loss term and atomic decay respectively. They are dimen¬ 
sionless and scaled with the typical resonance frequency is 
n = 51 GHz [^ . 7 is dimensionless, W is dimensionless and 
scaled with Q. 


We choose injection time 1/r = 1 /(Wxk) where k is 
the cavity decay rate and Nf,x is the number of atoms 
kicking the cavity field in the photon lifetime. The time 
elapsed when cavity is empty tq can be related to the 
interaction time such that tq = NemT. Thus, we can 
write 1/r = T{l + Nem)- Here, Nem is a factor introduced 
to measure tq in terms of r so that Nem = {NexHT) — 1. 

We solve the master equation by numerical methods 
and compare the results with the developed theory. We 
use QuTip package [S^ in Python software to solve the 


Figure 5. (Color online) Extracted work {W) and efficiency 
[rf) of optical resonator photonic Carnot engine, with A^ + 1 
level atom phaseonium (NLAP) fuel, depending on the num¬ 
ber of degenerate coherent ground state levels N for differ¬ 
ent decoherence factors (a)-(b) ^ = exp(— 2 ;), (c)-(d) ^ = 
exp(— A'a;), (e)-(f) ^ = exp(— A'^x), where x = 7 ^/ 7 . Coher¬ 
ence parameter is A = 10“® and the initial thermal coherent 
atomic temperature is T;i = 4 in units of x values 

are 0.15, 0.1, 0.05, 0.001 for (a)-(b), 0.14, 0.12, 0.1, 0.08 for 
(c)-(d) and 0.01, 0.008, 0.006, 0.004 for (e)-(f) in decreasing 
order for the lower to upper curves respectively. The param¬ 
eters g = 6.28 X 10“^, r = 8 X 10“®, k = 2.86 x 10“^ and 
7 = 4.68 X 10~® are the coupling frequency to the cavity field, 
atomic injection rate, cavity loss term and atomic decay re¬ 
spectively. They are dimensionless and scaled with the typical 
resonance frequency is = 350 THz [ 3 ^. rj is dimensionless, 
W is dimensionless and scaled with S2. 

master equation. We perform single atomic injection in 
two steps. First step is the atom-cavity field interaction 
(during r) and the second one is the free cavity field evo¬ 
lution (during tq). The master equation for the first step 
is written under Markov and Born-Markov approxima¬ 
tions as [ 3 ^ 

Af-l-l N 

p = -z[H,p]+7E^[^™] +(13) 

m n 

where last two terms stand foi^ure spontaneous emis¬ 
sion and pure dephasing [40l443l| . respectively. Here, 
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C[x\ = {2xpx‘^ —xx‘^p—px‘^x')l2 is a Liouvillian superoper¬ 
ator in Lindblad form and = \r){am\^ Lf^ = \bn){bn\- 
We include an auxiliary state |r) to the atomic state 
space to model the decay of the excited and the lower 
levels. Presence of |r) is not altering the initial phaseo- 
nium state. The auxiliary state is unpopulated and at 
the lower level energy. Its use allows for faithful simula¬ 
tion of the excited state and the degenerate ground state 
{am = a, &i,..,& 7 v) decay equations in Eq. (IA43I) . This 
decay model is already used in the original master equa¬ 
tion developed for the two-level phaseonium engine Q. 
Different decay models, for example decay of excited level 
to the lower levels are employed for other systems such 
as many atom superradiant Otto engine [20j and simi¬ 
lar effect of beating decoherence with scaling up coher¬ 
ence is found. Present contribution discusses the original 
photo-Carnot engine Q as well as the objection to its fea¬ 
sibility due to dephasing, phenomenologically described 
by factor ^ [2lj. Our introduction of 7 is an additional 
decoherence channel not included in Ref. [2l|. We have 
found that 7 ^ can be analytically expressed in terms of 
^ (see Appendix for details). The microscopical mas¬ 
ter equation approach describes both the original photo- 
Carnot engine coarse-grained master equations with 
and without @ dephasing and generalizes them to the 
multilevel case; as our mesoscopic master equation, Eq.[21 
does analytically. 

In the second step, cavity decay {k) is present during 
the time interval tq in accordance with the key assump¬ 
tions of micromaser theory and the corresponding 
master equation is 

p= -i[H,pf] + K.C[a]. (14) 

In Fig. [51 we present the thermalisation process of the 
cavity field depending on different Nex values by depict¬ 
ing the photon number versus scaled time. Physical pa¬ 
rameters are given in the figure caption consistent 
with the Rydberg atoms in a superconducting Febry- 
Perot cavity [s^- For low values of N^x which corre¬ 
sponds to large iVem, we have zigzag like curves and for 
high values of Nex, we have smoother lines. Average 
photon number converges to the theoretical value for 
Nem = 12x10^. Thermalisation time is much longer than 
the convergence rate of coarse grained master equation 
and the microscopic exact method is much more costly 
numerically. 

In Fig.|7i;a)-(b), we express the effect of the number of 
degenerate atomic ground state levels N of the coherent 
atoms against the dephasing and the cavity loss mecha¬ 
nisms. Horizontal dotted and dotted-dashed lines are the 
analitical values of average photon numbers ( n^, n) of 
the no loss case for each N. In Fig. |7Ka) when decoher¬ 
ence channels are open, average photon number saturates 
below the analytical values for iV = 2 in accordance with 
the argument that 2LAP phaseonium can not beat de¬ 
coherence [ 2 l|. Average photon number exceeds n for 
N = 4 in Fig. jT^b) by keeping all the other parameters 
same. Thus, we show that the decoherence can be beaten 
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Figure 6 . (Color online) Time evolution of average photon 
number during thermalisation process of PCE held while reg¬ 
ular injection of 2LAP depending on different Nex parame¬ 
ters. Nex values are 4500, 1500, 150, 100 and 50 in decreasing 
order for the upper to lower curves respectively. Coherence 
parameter is A = 10“®, initial held temperature Tf = 1 and 
temperature of thermal coherent atoms is Th = 2 in units of 
M}/kB- The resonant held frequency is D = 51 GHz, cavity 
quality factor is Q = 2 x 10 ^°, atom cavity held interaction 
time is T = 10 ps, atom decay rate is 7 = 33.3 Hz, atom 
dephasing rate is 7 ,/, = 3.3 Hz and atom cavity held coupling 
is g = 50 kHz. Time is dimensionless and scaled with Q. 


by using higher N. To make the effect more visible in the 
figures, we take larger coherence magnitude, A = 10“^. 

We also compare the consistency of effective field tem¬ 
perature Teff, harvested work (IF) and efficiency ( 7 ) ver¬ 
sus N in Fig. [5] (a), (b), (c) between developed analytical 
and numerical results when dissipation channels open. 
We observe a good consistency between numerical and 
theoretical results in steady state. 


C. Preparation of the NLAP 
and its energy cost 

Typical methods to generate quantum superposition 
states, such as pulse area, adiabatic passage, or STI- 
RAP techniques [^, , utilize optical pulses interacting 

with the atomic system to transfer an initial atomic state 
to a target one. The initial and target quantum states 
are known and hence one can easily determine the re¬ 
quired unitary transformation between them. Physical 
implementation of the required propagator is however a 
much more challenging problem than the calculation of 
the transformation matrix. An efficient strategy to sythe- 
size the transformation matrix is to decompose it into a 
product of matrices, representing interacting steps that 
can be implemented by using optical pulses coupled to 
the atom. 

An arbitrary N dimensional unitary matrix U(7V) can 










0 200 400 600 800 

nt 

Figure 7. (Color online) Comparison of time evolution of aver¬ 
age photon number in presence of thermal and coherent atom 
injection with N — 2 and N = 4. The horizontal dashed 
and dashed-dotted lines stands for the analytical and h 
values in absence of loss mechanisms. Solid and dotted lines 
stands for time evolution of and h respectively in presence 
of dissipation channels. = 4000 for both subplots corre¬ 
sponding to To = 90 fis. Insets magnifies the lines between 
Qt = 600 and Sit = 800. All the remaining parameters are the 
same with that of Fig. [S] Time is dimensionless and scaled 
with fl. 


be decomposed into so called N generalized quantum 
Householder reflection (QHR) matrices or — 1 stan¬ 
dard QHRs and a phase gate [H, [2^. A generalized 
QHR is defined by 

M(:z;<^)=I + (e*'^-l)|t/)(^| (15) 

where I is the identity operator and the |iz) is the nor¬ 
malized column vector with dimension N, same with the 
number of the pulses, and (f) is an arbitrary phase fac¬ 
tor. The decomposition of U(iV) in terms of generalized 
QHRs can be written as 

U = M(izi; (^i)M(^2 ; (j)N). (16) 



Figure 8. (Color online) Comparison of numerical results 
(dashed line) with developed theory (solid line) of effective 
field temperature Tgg, harvested work W and efficiency rj 
respectively depending on the number of degenerate ground 
state levels N. = 12 x 10® and corresponding tq = 2.01 
ps. All the remaining parameters are the same that of Fig. (6] 
Tgfj, W and rj is dimensionless and scaled with fl. 


For (j) = TT, Eq. (USD reduces to 'M. = 1 — 2\v){v\ which is 
the standard QHR. The interaction represented by each 
Householder matrix can be described by a propagator 
which can be determined by the Morris-Shore transfor¬ 
mation p^ . 

Our N-l-1 level atom coupled to N optical pulses in a fan 
shaped transition scheme, or so called N-pod model, is a 
generic model that is used to discuss generation of arbi¬ 
trary multilevel superposition states. Under the Morris- 
Shore transformation, the lower levels of the atom are 
grouped into a single bright level coupled to an effective 
single pulse and N-1 dark levels uncoupled from the op¬ 
tical pulses. The propagator is then easily determined 
in this Morris-Shore basis. Back transformation to the 
original basis gives the full propagator, or the general¬ 
ized QHR matrix. Both the number of QHR steps and 
the number of pulses used in each step are in the order of 
N, therefore the total number of pulses to be used to gen¬ 
erate the target state would be in the order of N'^. This 
shows that the energetic cost of preparation of the target 
state scales with same with the work and efficiency 
scaling in the corresponding photonic Carnot engine. 

The preceding discussion is applicable to the case of 
mixed states as well, for which the normalized vectors of 
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generalized QHRs are defined as [28j 
1 / 2 sin ((/>/ 2 ) 


W^) = 


- 1 


I 1 ^7/ 


|w,)-|ei)). (17) 


Here, m is the column of U(iV), |ei) = [0,..0]'^ 
and (j) is an arbitrary phase where (pi = 2 arg(l — ua)— tt. 
It’s shown that for an N-Pod system any standard QHR 
M(z/) can be realized by single pulses with an rms pulse 
area H = 27r [2^. The corresponding unitary transfor¬ 
mation can only link the mixed states with the identi¬ 
cal dynamical invariants. In our case we consider initial 
thermal states out of a hohlarum transformed to a final 
state with small coherences. The initial and final states 
would then possess different spectral decompositions so 
that they cannot be unitarily connected. A resolution 
to this is suggested to exploit decoherence channels such 
as spontaneous emission or pure dephasin g, in combina¬ 
tion with the unitary transformation [2^, |^. We will 
not follow this route but use an alternative, which allows 
for a fully unitary procedure to generate desired coher¬ 
ences. As the coherences contribute additively, to exploit 
their scaling advantage we do not need an exact state but 
an approximate one would be sufficient. Accordingly, 
we can simply consider an approximate approach and do 
not specify an exact value for the coherences. We only 
need to keep them small enough to ensure slightly out of 
thermal equilibrium final state. We illustrate our strat¬ 
egy for = 2 case, and suggest that in principle larger 
NLAP can be generated by straightforward extension of 
this technique. Unitarity of our procedure also makes the 
details of generation process immaterial for the cost es¬ 
timation. The cost would be the same for other unitary 
equivalent processes to generate same states. 

Initial state of the atom out of the hohlarum at Th is 
the thermal density matrix 


Pc ~ UpthU^. The unitary transformation U can be de¬ 
termined from the matrix that diagonalize pc and found 
to be 

/I 0 0 \ 

U = 0 -0.707 0.707 . (20) 

\ 0 0.707 0.707 ) 

Writing pc = UpthU^, the fidelity between pc and pc is 
determined by T{pc,Pc) = \Tr and found 
be J^{pci Pc) — 1- Thus, U can be approximately be used 
as the unitary operator linking the initial and the target 
density matrices. 

U can be synthesized by using two standard QHRs and 
a phase gate as 

U = M(i.i;(^i)M(z/2;</>2)M(:/3;</>3) ( 21 ) 

wherec))! = p 2 = p 3 = and M(r' 3 ; ^ 3 ) = $( 0 , 0 , $ 3 ) is a 
one-dimensional phase gate. By using Eq. (II3 one finds 
the normalized column vectors to be \vi) = [0,0,0]^, 
\v 2 ) = [0,0.924,-0.383]^ . The example for = 2 
here illustrates the basic principles to generate arbitrary 
NLAP. One can use N pulses for each N — 1 standard 
QHRs and a phase gate to synthesize a unitary trans¬ 
formation matrix, which is the diagonalization matrix of 
a quasi equilibrium thermal state at Th with small co¬ 
herences. When the pulses applied to the actual atom 
out of hohlarum at Th-, the final state will be approxi¬ 
mately the same with the target state used to determine 
the properties of the pulses. 

We can now estimate the energy cost Uc of generating 
Pc using Uc = N'^Up where N'^ is the total number of 
pulses used in the QHR technique and Up is the energy 
of a single pulse. Up can be determined from the pulse 
area A as in [ 2 ^. For a square pulse of duration Tp 
and amplitude Ep we have A = dEpr/h, where d is the 
magnitude of the dipole moment 


, w-i-i 

Pth=^e-^^= (18) 

where /3 = l/fcsT/i inverse temperature [ks = 1 ) 
and Z = Tre“^^ is the partition function, with 
H = being the atomic Hamiltonian. 

Taking N = 2 and Th = 2, we hnd pth = 
diag(0.327,0.384,0.384). Target density matrix pc is 
taken to be 

/ 0.327 0 0 \ 

Pc=\ 0 0.384 0.000001 , (19) 

\ 0 0.000001 0.384 j 


where e is the vacuum permittivity and c is the speed of 
light. Taking the pulse area A = 27r we hnd 

2TTh I , 

The intensity of the pulse is given by Ip = c€o\Ep\‘^/2. 
The pulse energy in a beam of radius r;, can be estimated 
by Up = nr'^IpTp. Using H = 2'KcfX and p = Xj2Trrt, 
where A and ( are the wavelength of the optical held and 
the radial beam divergence, respectively, we hnd 


where the off-diagonal elements between degenerate 
ground state levels are taken real and much smaller than 
the diagonal elements that are equal to those of the pth- 
The initial and hnal density matrices have distinct dy¬ 
namical invariants in their spectral decompositions and 
hence they cannot be linked by a coherent evolution. 
Let us assume however an approximate link such that 


TT^ 1 1 

Up = hn——-. 

6 Tp-fC 


(24) 


Taking 1/Tp^ ^ 2 and C ^ 0.5 [l^, we hnd Up ^ 12hQ. 

Total energy cost to reach steady state per cycle is 
Uss = mUc = mN'^Up where m = rAtg. Here, m is 
the number of atoms needed for thermalisation, r is the 
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injection rate, and At^ is the time elapsed to reach the 
steady state. In our results we have found that the har¬ 
vested work per cycle is much less than the resonance 
energy, W <C fill, thus the generation cost of NLAP fuel 
is several orders of magnitude larger than the harvested 
work Uaa ^ W, which confirms that the second law of 
thermodynamics obeyed in photonic Carnot engine with 
NLAP. Generation cost is not included in the thermody¬ 
namic efficiency but it can be significant figure of merit in 
the round trip efficiency. To make such photonic Carnot 
engines more appealing for certain applications, it is nec¬ 
essary to increase their round trip efficiency as well. For 
that aim one could consider the cases of larger compres¬ 
sion ratios (Afl ^ fl) and high operation temperatures 
[ksTh ^ Ml). Our focus here is on the discussion if such 
engines can operate under decoherence. Despite the neg¬ 
ative conclusions for two level phaseonium case [2lj, we 
have found that larger phaseonium fuel allows for opera¬ 
tional photonic Carnot engines. The question of how to 
increase their round trip efficiency requires further anal¬ 
ysis which is beyond the scope of present contribution. 


IV. CONCLUSIONS 

Summarizing, we examined scaling of work and effi¬ 
ciency of a quantum heat engine with the number of 
quantum resources. Specifically, we considered a pho¬ 
tonic Carnot engine with a multilevel phaseonium quan¬ 
tum fuel. We derived a generalized master equation for 
the cavity photons, which forms the working fluid of the 
engine, and determined the steady state photon number 
to calculate the work output and thermodynamic effi¬ 
ciency. We find that they scale quadratically with the 
number of quantum coherent levels N. We examined 
the case of degenerate levels to get analytical results and 
to examine scaling laws against decoherence due to cav¬ 
ity dissipation and atomic dephasing. We verified our 
analytical results with detailed numerical methods and 
have shown consistency of coarse grained analytical re¬ 
sults with the microscopical numerical approach. Gen¬ 
eration of multilevel phaseonium fuel using Morris-Shore 
transformation determined quantum Householder reflec¬ 
tion technique as well as its cost are examined. Using 
typical parameters in modern resonator systems, such 
as circuit QED, our calculations reveal that decoherence 
due to cavity dissipation could be overcome by the mul¬ 
tilevel quantum coherence even in the presence of large 
dephasing rate. If the dephasing rate increases with V, 
then work and efficiency can still overcome the decoher¬ 
ence and retain their scaling up to a critical number 
of coherent levels. 
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APPENDIX 

We generalize micromaser mesoscopic master equation 
treatment @ applied for three-level phaseonium engine 
U to multi-level case. In ALAP model, the Hamiltonian 
of the whole system is H = Ho + Hf where label k implies 
injected k*^ arbitrary atom. We adopt the notation of ref. 
[44| for the Hamiltonian and relevant quantities. Here, 
Hq = hua\a'"){a’^\ + and 

N 

= hg^\a^){b\\ae-'^^^ + H.c. (A25) 

i=l 

in the interaction picture where wi, are atomic energy 
levels, g is the atom-field coupling coefficient and D is the 
single mode cavity frequency. Here, we assume all levels 
coupled to the excited one with the same coefficient g. 
The equation of motion of overall system is 

p =-^[H,p\+C a[p]+C f[p\ (A26) 

where La [p] and L / [p] are the Liouvillian superoperators 
expressed in the main text corresponding to atomic and 
field degrees of freedom respectively. 

The equation of motion of the radiation field which is 
the working substance of the heat engine can be found 
by tracing out atomic part as 

Pnn = -;^^(Trat[i^^p'=]„„) 
k 

+Lva.tLA[p]nn + Trat/:/[p]„„ (A27) 

where pnn = {n\p\n). Here, [H,p\nn = {n\{Hp - 
pH)\n) = T.ni'<^\H\m){m\p\n) - {m\H\n) = 

J2m.i^rimPmn- PnmHmn)- In micromaser theory, due to 
the short atom-cavity interaction time, last term, cavity 
decay is usually ignored when the atom is inside the cav¬ 
ity. The second term is treated perturbatively and will 
be considered for zeroth order in g. Here, we will keep 
it but assume it can be treated independently. One may 
write the partial trace operation over atomic degrees of 
freedom for a random single atom as, 

Tia.t[H,p]„n = '^{a,n\[H, p]\a,n) 

06 

N 

= {an\[H, p]\an) + ^{b^n\[H, p]\b,n) 

(A28) 

where a are the atomic basis as expressed at the right 
hand side. Each term of Eq. (IA28I) can be calculated 
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by using the selective rules of the Hamiltonian between 
certain levels n and m; for instance inserting for the first 
term of Eq. (IA28I1 we have, 


{an\[H, p]\an) 


(an\{Hp — pH)\an) 
'y^{{an\H\a'm) {a' m\p\an) 

<y'm 

{an\p\a m){a rn\H\an)'\ 




ct'm^an Pan,a'mH q^' rn,an\■ 

(A29) 


We use pf = p hereafter for simplicity. Here, 

Cf[p]nn = (n|f (2apa'l' - pa'^d - d'^dp)\n) 

= K{{n + l)pn+i,n+i-np„n}. (A35) 

In order to proceed the calculation, any single term in 
the summation of Eq. (IA34I) should be calculated and in¬ 
serted therein. The terms can be obtained by the integra¬ 
tion of corresponding equation of motions by using the 
selecetive rules of the Hamiltonian as expressed above. 
We evaluate atomic equation of motion 0*^ order in g 
first, we have 


Hi is the hamiltonian Eq. (IA25I) to be inserted into 
Eq. (jA29|l which can be written as. 

Hi = ngR+de-^^*+ hgR_d^e^^\ (A30) 

where R+ = R- = Expressing the 

terms of Eq. (jA29ll conveniently, we write 

{R+)aa> = (a|(|o)(6i| -I- ... -I- |a)(6Ar|)|a') 

= (Sb^a' + + dbj^a') (A31) 

and {R-)aa> = 0. Besides, dnm = {n\d\m) = ^/m{n\\m — 
1) = Substituting these terms into the first 

part of Eq. (IA29I) we have Y.a'm{Ha.n ^OL'mPoi'm,an\ — 
hge~'’^*-y/n + l{p b^n+i,an + ■■■ + PbMn+i,an)- The second 
part of Eq. (IA29I) is simply the complex conjugate. Hence 
the hrst term of Eq. (IA28I) is 


N 

{an\[H,p]\an) = hge~"^W n -H Pbin+i,an - c.c. 

i=l 

(A32) 

The second term of Eq. (|A28I) would be cal¬ 
culated by similar considerations. We can write 
{bin\[H, p]\bin) = In this case 

{R+^biot' — {,R—)bia' — (n||oi) — and d^^ — 

(n|a’^|m) = ^/rnH^{n\\m -I- 1) = ^/m+ lSn,m+i- Then, 
{b^n\[H, p]\bin) = hgef^^ynpan-iMn-C-c. and hnally the 
second term of Eq. (IA28I) becomes, 

N N 

'^{bin\[H, p]\bin) = hge''^*'y/n'^ Pan-iMn - c.c. 

i=l i=l 

(A33) 

Inserting these results into Eq. (IA27|) and after some re¬ 
arrangements we have the field equation of motion, 

N 

Pun = -9 + Ie“*^* ^ pl,n+l,an 

k 

N 

- iy/ne-^^* ptn,an-l) + c-c.} 

(A34) 


Pa — ^ PaI ^ ^ 

cx,oc' 

+ (A36) 

i 

where La^a' = |ck)(Q^^I and Lb^^bi = \bi){bi\. Final two 
terms of Eq. (|A36p corresponds to Ca [pa] with a' ^ a 
and i = l,..,iV. Here, a = {a, bi,bmj and ja,a' is 
taken equal to 7 for simplicity. Note that the atomic 
part of the master equation is for the case of pure de¬ 
phasing & relaxation and we follow the usual assumption 
of micromaser theory [ 2 ^ that cavity decay and atomic 
dynamics can be separately treated. When atom is inside 
the cavity, decay is not included. 

The equation of motion of the term of the first 
summation of Eq. (IA34I) for a single atom is 

Pb^n-t-l.an — ig^TL -\- Ic {pan,an 

{Pbin+l,bin+l “t“ .. 4- Pbin+l,b]\^n+l')}• (A37) 

Here, we have neglected the the matrix element {bin + 
l|£/[p]|an) in accord with the assumptions indicated 
above. The equation of motion for pb^n.an-i which is 
the A term of second summation of (IA34I) could be ob¬ 
tained by simply replacing n —>■ n — 1 in Eq. (IA37I) . that 
is 

PbiU^an—l — igy/nC {pan—l.an—1 

(^Pbin,bin 4“ .. 4- PbiJi^b^n')} (A38) 

Therefore we obtain Pbin+i,an Pbin,an-i terms by in¬ 
tegrating Eq.s ljA37l and IA38P formally in the following 
form. 


PbiU+l 


^^ = -igAr+T [ 


^{Pan,an {Pb in+l,bin+l 4- .. 4- Pbin+l,bNn+l')^ 


(A39) 


PbA 


,au-i = -wV^ f 

Jtkr, 


\Pan—l,an—l {PbinAin 4- 4- Pb 

n,b\Tn) } 


(A40) 


+ ^f[Pf]nn- 
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The terms of Eq.s (IA39I and IA40I) can be factorized to in which the solutions are 
atomic and field density matrices for zeroth order solu¬ 
tion in g. For instance, 


Pcm,ani^ ~ Pa,ai^ itko)Pn,n{t ), 

Pan+l,an+li^ j ^feo ) ~ Pa,ai^ itko)Pn+l,n+lit )) 


Pb°n,bjni^ J^kg) — Pb°,bji^ , tkg)Pn,n{t ), 

Pbin+l,bjn+A^ i^fco) = Pb°,bj^^ : ^fco )Pn-|-l.ra+l ) 

(A42) 

Here, which are the initial atomic density ma¬ 

trix elements, obey the respective atomic equations of 
motion 

Pi°ait',tkg) = 

Pb°,bj tko) = -(iWb.b- -h 7 + l<j>)pt%j (i': tkg) 

(A43) 


°^ptlbAtko,tko) (A44) 


„ko _ „-(i<^i>ii>,+7)(i-tfco) „fco 
Pbi,bj — ^ Pbi, ^ 

where 7 = 7 -I- 70 . Eq.s (IA43I) and (IA44I) imply that 
excited and ground state levels decay to a lower level. 
(A41) off-diagonal elements of the atomic density matrix 

are equal to (4o,4o) = The is 

assigned with the coherence preparation. Thus, we can 
’ express Eq.s (IA39l and lA40ll by using zeroth order atomic 

Eq.s (|A42IIA44ll to find first order solutions in g. 


pL+i,an = -wV^^ f 

Jtkg 

f) _ p 

HaaHnn ^ Hb 


X{e ‘'“oVaaPnn-e ‘"o) p^ibi Pn+ln-Hl 




(A45) 


Likewise, 


Pb. 


JUr. 


-7(*'-*fco)^0 

HbibiHnn 




'tk, 


|Pbi,bJe*'P*V„„}. 


(A46) 


Putting all these results into Eq. (IA34I) we have the field equation of motion, 

Pnn = - 9 ^ E f + 1) f; ,,p„+i„+i 

k *^^^0 

N 

Vb.,bJe'‘^'"Pnn)] +C.C.} 


i=l 

_ g(.i‘^bibj+l){t'-tko)\^0 ^ |pZ0 


(A47) 


Before proceeding, we replace the summation over number of injected atoms by integration over injection time as 
^/_oo *^^0 where r is the injection rate. We also define = uiabi — where uiabi = bOa — uibi- Then, 

N 


Pnn = -rg‘^ [ dtko [ -k l)p„n - np„_i„_i) 

d-oo Jtkg 


N 


- ,J((n + l)p„+i„+i - np„„) 


i<j 

X ((n -I- I)p„+ 1„+1 - npnn) + C.cj -I- Cf[p\nn- (A48) 
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Note that = e ^ *'=0) |e *'^0 while pb^&j = e +T'h* **“0)|p0^ |e* 0 o, Evaluating the integrals 

in (IA48I1 over t' and tka after changing integration order as /* dtka C dt' = /* dt' /* dik^ we have 

Kq 


(.N N 

Pnn = -rgH -(E _ia.. + -,. )'’°d(" + iK~- 


7 -iAi + 7 

i:((y 




i<j 


{-iAi + 7 ) (iuJbibj + 7) (-*Aj + 7 ) {-iujbibj + 7) 


7 -iA^ + 7 / 


X [{n + l)pn+ln+l - npnn] + C.C. > + Cf[p]r 


(A49) 


We proceed by summing each term with their respective complex conjugates, some arrangements then, 


r 1 Af 27 1 ^ 

Pnn = -rg^ “I “ ( X/ /\2 ^2 ) Paa [(^ + ^)Pnn - «Pn-ln-l] ” “ ( X/ 


-jO 

^ ^2 JPbib^ 


7 V^ A, +7 


+ 


E[( 


2 cos 4>ij{AiOJbibj + 77) + 2 sin 4>ij{u^b^bj - Ai^) ^ 2 cos (jiij (77 - AjUb^bj ) + 2 sin (j)ij{uibibj + Aj^) ^ ^ ^ 


i<j 


(A|+7^)K-,,+7^) 




X [{n + l)p„+i„+i - npnn] 


\Pbih 
(A50) 


and by using Eq. (IA35I) . finally we have, 

Pnn — [^(n A l)pnn ^Pn — In — l] A (^Rqq A d^gc^ ^ \jbPnn A l)pn+ln+l] ^ 

A K{(n A l)p„+l,n+l - »^Pnn} (A51) 


where 


N N S 

i=l ® ' i=l i<j * 


(A52) 


2 cos (/)ij {AiWbibj A 77 ) A 2 sin {ujb.bj - Ai^) 2 cos (l)ij (77 - AjUJbibj ) A 2 sinA Aj 7 ) 




(A|A72)(A 


bibn 


■7^) 


(A?+7^)K,, 


+ 7^) 


(A53) 


Aij- = Wabi ,3 —oJabij = bJa — bOb^j and R = rg"^. Note that Rg^ has N number of terms and Rg^ has S = A(A— l)/2 
number of terms in the summation. Since we seek the solutions in the steady state, we obtain the steady state photon 
number by solving = 0 where n = J2n ’’^-Pnn and we write by using previously obtained pnn as 


h = -RKa Paa E Ti{tI + RKaPaa ^ ^ Pn—l,n—l ^ ^ ^ Pnn ^ ^ Tl(^Tl + 

n n n n 

-rr^pT. Pnn H“ rr^pT. n{n A l)pn+i,n+i A k ^ ^ n{n A l)p„+i^n+i ~ k, ^ ^ Pnn- (A54) 


Then, we insert n —>■ n — 1 for and n —> n A 1 for 
Pn+in+i and pn-in-i terms respectively so that we get. 


n — RKaPaa ^ ^ (n A l)pnn RRqq ^ ^ ^Pnn 

n n 

— RRg^ ^ ^ Pnn 
n 

RRaPaa{j^(fi “t“ 1) R^tpi^RgQ ^gc) blTlfp (A55) 
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Solving = 0, we have 


_ _ ^aPaa 

^ -Rgo + ~ KaPaa 

1 

^80 I Rqr I_ K, 

KaPaa KaPaa RKaPaa 


^liRgol^aPaa — 1 ) which is the average photon number 
in the absence of atomic coherence and in the absence of 
cavity decay k. 

Finally we look at the degenerate ground state case 
i^abi^j = 0, Aij = 0). In this case, h can be simplified 
to h = Pe/{Pg - Pe) where Pg = Paa, Pg = Pbib, for any 
i. The simplihed forms of other parameters are 


We write the final form of the steady state photon num¬ 
ber after some arrangements, 


where 


= 


1 + n. 


R„ 


KaPaa 


(A57) 


Kg 


^2 ’ -^90 


2NP, 


Rpc 


2N{N -1) cos (pX 
77 

(A59) 


for 9 = TT. Analytical decoherence term for degenerate 
case can be identified in Rg^ expression such that 




1 -I- n 


RKaPaa 


(A58) 


Here, is the average photon number in the absence of 
atomic coherence in terms of average photon number h = 


for 70 < 7 . 


e = (i + ^)-i 
7 


^ g-7.#./7 


(A60) 
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